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We propose a microscopic physical mechanism that stabilizes coexistence of the Kondo effect and 
antiferromagnetism in heavy- fermion systems. We consider a two-dimensional quantum Kondo- 
Heisenberg lattice model and show that long-range electron hopping leads to a robust antiferromag- 
netic Kondo state. By using a modified slave-boson mean-field approach we analyze the stability 
of the heavy antiferromagnetic phase across a range or parameters, and discuss transitions between 
different phases. We also address connection to experiments on heavy fermion compounds. 
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Introduction. The study of complex phenomena ex- 
hibited by materials with competing ground states is the 
primary focus of the modern condensed matter physics. 
In systems where local magnetic moments interact with a 
conduction band there are two opposing quantum many- 
body effects: Kondo screening (formation of singlets be- 
tween the local moments and itinerant electrons), and a 
long-range magnetic [often antiferromagnetic (AF)] or- 
der. Since metallic Kondo phase and magnetic ordering 
involve the same local-moment degrees of freedom, they 
contest the same entropy. This competition is at the 
heart of the rich variety of phases observed in heavy- 
fermion (HF) /-electron materials under tuning of exter- 
nal parameters like pressure, doping, or magnetic field 

In the Kondo screened phase the Fermi surface (FS) 
volume accounts for the local spins. Quantum (zero- 
temperature) phase transitions between this HF metal 
with large FS and small-FS magnetic states remain an 
actively debated subject j2|!3]. Of particular interest is 
whether the Kondo screening collapses precisely at the 
onset of magnetism [H [S] or the magnetic transition is of 
a spin-density wave (SDW) type with no concomi- 

tant change of the FS volume. Both scenarios are found 
in experiments [3]. In Ce3Pd2oSi6 [8] Kondo screening 
disappears inside a magnetically ordered phase, while in 
CeCu2Si2 the AF transition is likely of the SDW type 
[S]. In YbRh2Si2 data indicate that Kondo screening 
disappears at the magnetic transition IIQI [TTj , while un- 
der Co and Ir doping the two transitions separate 
Developing microscopic theories which exhibit such mul- 
titude of phases proved difficult pQ |31 [H] . Recently it 
was conjectured [31 [13] that frustrated magnetic inter- 
actions between local moments give rise to a variety of 
phase that include magnetically ordered as well as para- 
magnetic states with both large and small Fermi surfaces, 
but detailed theories are still lacking. 

In this Letter we propose a microscopic mechanism 
which controls the coexistence of the Kondo effect and 
AF long-range order. We consider the two-dimensional 
Kondo-Heisenberg lattice model with short-range AF 
interactions between local moments, and conduction- 



electron hopping beyond nearest neighbors (NN). By em- 
ploying the notion of a "spin-selective Kondo insulator" 
[141 115) introduced in the context of ferromagnetism in 
Kondo lattices we show that in a bipartite Kondo lat- 
tice with NN and next-nearest neighbor (NNN) electron 
hopping the same physics leads to a robust AF Kondo 
["K+AF" in Fig. [^a)] phase. The stability of this state 
is controlled by the relative magnitude of short- and long- 
range hopping amplitudes. The AF Kondo phase has a 
large FS, and is separated by a second-order quantum 
phase transition from the small-FS AF metal. In con- 
trast to previous studies [TH [THl [I7| we focus on the ex- 
perimentally relevant regime away from half-filling, and 
perform an unbiased calculation for competing orders. 

Model and approximations. The essential physics of 
magnetic HF materials is contained in the Kondo- 
Heisenberg lattice model (KHLM), 

H =-ffcond + -ffKondo + ^^Hcis = (1) 

{ij),a i,a 

+ Jk ^ SiSi + Jh ^ SiSj. 
I {ij) 

This Hamiltonian describes a system of itinerant elec- 
trons, Cia with spin indices a, j3 =ti interacting with 
local spin- 1/2 moments Si via the AF Kondo coupling 
Jk- We consider this Hamiltonian on a square lattice 
with sites labeled by i and j. The competition to the 
Kondo screening is provided by the AF exchange J^j > 
between NN spins. The hopping amplitude tij — t > 
for NN hnks {ij) = (ij), and tij = t' for NNN sites, 
[ij) — {{ij)); for all other sites tij = 0. The electron spin 
operator is S; — (7apcl^Cii3/2, where or are the Pauli ma- 
trices. The chemical potential fic controls the conduction 
band filling. 

We employ the hybridization mean-field (HMF) ap- 
proach [T], where local spins are represented in terms of 
spin-1/2 pseudo-fermions. Si — (Tapfl^Sip/^, subject to 
the constraint flafia — Ij which removes unphysical 
empty and doubly-occupied states from the single-site 
Hilbert space. The core idea of HMF is to treat this con- 
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straint on the average by introducing the pseudo-fermion 
"chemical potential" fij. The interactions in Eq. ([!]) 
now contain four fermion operators and are decoupled 
within the Hartree-Fock approximation. We consider de- 
couplings in all possible channels, namely: 
(i) Magnetic channel 



(a) 



^^Kondo + ^^Hcis 



(2) 



i (ij) 

where magnetic order parameters (OPs) are defined as 
Mi = (Si) and rrii — (s^), and the arrow indicates that 
we omitted c-number terms. 

(ii) Pseudo-fermion dispersion ( "spin-liquid" ) channel 

^Hcis ^ XI '^"'"'^/J'/3 [(/L'/i/J>/jV/»a+h.c.] . (3) 



(iii) Kondo hybridization channel 



— 2^ XloX^O +^2^Xl Xi 

i 



(4) 
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with = '^afifia^^p/^'^^ A* = 0...3 and cr°^ = Sap. 
Physically, the x-operators are Schwinger bosons [1"51 
which create local singlet and triplet states resulting from 
the Kondo coupling between localized and itinerant spins. 
In the usual picture of Kondo singlet formation only xo~ 
boson is present. The magnetic order admixes other com- 
ponents, so that the x-representation of Eq. Q captures 
the SU{2) invariance of the Kondo interaction [T9". 

Since we consider only uniform and commensurate AF 
states on the square lattice, there are two wavevectors 
in the problem: q — Q and q — Qo = (7r,7r). In the 
AF phase all OPs are site-dependent. Thus the pseudo- 
fermion density n( = ^aflafia will acquire an un- 
physical spatial dependence. To suppress the Qq har- 
monic Uq^ we impose an additional constraint, H — >■ 

H—fJ-Q J2i e^^°^^n{, where fiq is the Lagrange multiplier. 

In the rest of the paper we study the phase diagram of 
the KHLM as a function of t' , Jh and temperature, T, at 
fixed density n'^ = 0.8. We choose the spin quantization 
axis along the z-direction and omit the corresponding 
vector indices whenever possible. Then the HMF equa- 
tions are solved numerically on the 32 x 32 square lattice 
with periodic boundary conditions. 

Origin of the AF Kondo phase: t — Q limit. It is 
known [TU [TS] that the Kondo lattice Hamiltonian with 
Jh = and NN hopping exhibits a ferromagnetic HF 
state characterized by finite Kondo hybridization and 
magnetic polarization. In that "spin-selective Kondo in- 
sulator" phase IHj, away from half- filling, all minority- 
spin and the corresponding fraction of the majority-spin 
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FIG. 1. KHLM with t = Q. (a) T = phase diagram: K 
- singlet Kondo phase, M - local moment magnetic phase, 
MKi,2 - magnetic Kondo phases. M and MKi,2 are ordered 
at q = if Jh < and (tt, tt) for Jh > 0. All phase transitions 
are 1st order. Inset: Square lattice with t = 0. The arrows 
denote local moments, (b) and (c) Hybridization (xo.z), mag- 
netic (A/ and m) OPs, and total magnetization M + m for 
t'/Jx = 2/15. The plateau in (c) reflects Eq. 



conduction electrons screen part of the /-electron spin, 
while the magnetic moments of the remaining conduction 
and local /-electrons are antiparallcl to take advantage of 
the Kondo coupling. To reveal the physical mechanism 
responsible for coexistence of the Kondo effect and anti- 
ferromagnetism, it is instructive to start with a limiting 
case with NN hopping t = 0, but finite NNN hopping t' . 
The model Eq. ([T]) then reduces to two interpenetrating 
square Kondo (sub)lattices, coupled by the Heisenberg 
term [see inset in Fig. [ija)]. 

When Jh — and t' / Jk is large enough, each sub- 
lattice enters a ferromagnetic Kondo state with non-zero 
{Xo,z) and magnetic OPs. The ground state of this sys- 
tem is continuously degenerate with respect to the angle 
between the magnetizations of the two sublattices. Fi- 
nite Jh lifts this degeneracy and stabilizes long-range 
magnetic order coexisting with Kondo singlets across the 
entire lattice. Since for t = only the Heisenberg interac- 
tion couples the sublattices, the ground state energy is an 
even function of Jh, and FM (AF) states are stabilized 
for Jh < (Jh > 0). Below, m and M denote uniform 
(for Jh < 0) and staggered (for Jh > 0) z-axis magne- 
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tizations in the c- and /-channels respectively. The vec- 
tor part of the hybridization, (xz), has the same Fourier 
components as m and M, while the singlet hybridization 
amplitude (xo) is always uniform. 

Fig. (TJa) shows the T = mean field phase diagram of 
the KHLM with t = 0. For small t' and \Jh\ the system 
resides in the Kondo singlet state. In the opposite limit, 
the Kondo effect is suppressed in favor of magnetism. In 
the intermediate range of parameters several magnetic 
Kondo states (MKi 2) are stabilized. All phase transi- 
tions in Fig. [ija) are discontinuous as illustrated by the 
I Jij [-dependence of the hybridization and magnetic OPs 
in Fig. ^h). The states MK1.2 can be distinguished 
based on the commensurability condition |14j 



(a) 



2(M + to) 



(5) 



satisfied only inside MK2. Consequently, the state MK2 
has a plateau in the total magnetization, shown in Fig. 
[ijc). We also note that the phase MKi is quite fragile 
and may be an artifact of the HMF approximation. 

Coexistence of antiferromagnetism and Kondo effect. 
The above picture survives in the physically relevant limit 
^ i'/t ^ 1. Similar to the t = Q case, when Jh = and 
Jk are small compared to the electron bandwidth the 
system undergoes a transition to a ferromagnetic Kondo 
state. To avoid this, we fix the Kondo coupling at a fairly 
large value Jk = Qt, so that for = and any t' G [0, 1] 
the ground state is the non-magnetic uniform HF phase. 

Numerical solution of Eqs. ( 1 l-Q exhibits a complex 
phase diagram shown in Fig. 2 a). Its most important 
feature is the existence of a critical point D at t'^/t ~ 0.46. 
For t' < t'^ there is a direct 1st order transition between 
the large-FS HF and small-FS AF states. On the other 
hand, for t' > t'^ there is a large parameter range (shaded 
region in the figure) where Kondo screening coexists with 
antiferromagnetism. The hybridization OP is uniform 
in the K state, (xio) = (xo)j and vanishes in the AF 
state. The latter is characterized by = Me"5oa;i 
and mf — toc''^"^* . In the coexistence phase (K -I- AF) 
the singlet component of the hybridization acquires a 
small staggered component, {xm) = ixo) + (xo)e"^°^' 
(l(xo)/(xo)| ~ 0.01), while the vector part of the Kondo 
hybridization is purely staggered, (xiz) = (Xz)e"^°^S 
tracking the spatial variations of the magnetization. 

The coexistence phase is separated from the AF state 
by a 2nd order transition. However, its boundary with 
the pure Kondo state is more complex. In Fig. [2]^a) 
the line AB is the 2nd order transition, line DCB is 
weakly 1st order, and BD and BE are strongly 1st order 
transitions. The AF Kondo state inside regions ABE 
and BCD differs from the rest of the intermediate phase 
only in its value of the staggered magnetization M + m, 
see inset of Fig. ^a). Absence of the magnetization 
plateau in this phase [cf. Fig. [ijc)] indicates that the 
commensurability condition ([5| no longer holds in the 
presence of the NN hopping. 
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FIG. 2. KHLM at T = 0. (a) Phase diagram with K - 
singlet Kondo phase, and "K -I- AF" - large-FS AF Kondo 
state. The dashed line is a continuation of the boundary 
between K and AF phases, obtained if their coexistence is 
prohibited in the calculation. Thick (thin) lines denote 1st 
(2nd) order transitions. Inset: Uniform part of the Kondo 
hybridization OP and total staggered magnetization for t' /t = 
0.81. (b) Evolution of the heavy quasiparticle FS with t' for 
Jh = 0. The Lifshitz transition at t' /t ~ 0.5 corresponds to 
point C in panel (a). 



The kink in the phase boundary at point C is due 
to the change in the topology of the heavy quasiparti- 
cle FS (Lifshitz transition) at t' /t - 0.5 in Fig. ^h) [2D]. 
Hence, the intermediate phases ABE and BCD may sig- 
nify tendency towards incommensurate ordering, rather 
than the (tt, tt) AF state that we consider. More stud- 
ies are needed to determine the exact spatial structure 
of the coexistence state. Although the AF Kondo state 
exists only because of NNN hopping, increasing t' even- 
tually destroys it when Jk becomes small compared to 
the bandwidth. 

Finally, in Fig. [3] we present the quasiparticle density 
of states (DOS) in the three phases along the line t' /t = 
0.81 in Fig. [SJa). As expected, the DOS in the AF Kondo 
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FIG. 3. Quasiparticle DOS for a fixed t' /t = 0.81 and Jh 
inside the singlet Kondo phase (left), magnetic Kondo state 
(middle) and AF metal (right). The energy e is measured 
relative to the Fermi level. In the left and middle panels the 
Kondo peak is apparent. The two peaks in the right panel 
correspond to flat /-levels split by the exchange field in the 
magnetic phase. 



phase shows a peak near the Fermi surface, similar to the 
singlet Kondo state, and is very different from the DOS 
in the small-FS AF metal. 

Finite temperature behavior. To assess thermal sta- 
bility of the AF Kondo state, in Fig. |4f|a) we present the 
finite-T phase diagram of the KHLM computed along the 
J|-shaped path centered at point D of Fig. [2][a) . In the 
HMF analysis the Kondo screened phase (K) always van- 
ishes via a 2nd order phase transition to a small-FS metal 
at a critical temperature Tfi-( Jjj, i'). Depending on Jh 
and t' the HF metal may become unstable towards either 
a pure AF state or a "K -I- AF" phase [shaded regions in 
Fig. ^a.)] at a Neel temperature TN{JH,t') < T^. The 
latter phase extends over a noticeable temperature range 
with Tfq/Tc ~ 1/4 and are separated from the singlet 
Kondo state by a 2nd order transition. This is illustrated 
in Fig. |4]jb) which shows the temperature dependence of 
all OPs. In contrast, pure AF phase is separated from 
the HF state by a 1st order transition. 

Discussion. Our study identifies the long-range elec- 
tron hopping as a physical mechanism responsible for the 
robust coexistence of antiferromagnetism and Kondo ef- 
fect in the description of HF materials via the Kondo- 
Heisenberg lattice model. The parameter t' /t controls 
whether the Kondo screening can vanish precisely at the 
onset of the magnetic order, or via and intermediate coex- 
istence regime, covering the range of experimental data. 
Although we considered a square lattice, our results can 
be straightforwardly applied to any bipartite graph. Our 
phase diagram presented in Fig. [2|^a) can be verified ex- 
perimentally in HF systems by changing t' /t with pres- 
sure or chemical doping. 

We emphasize that our study is fundamentally differ- 
ent from previous works, notably Ref. [T^] which con- 
sidered the model ([T]) with t' = and found a phase 
diagram that included the heavy SDW state similar to 
our AF Kondo phase. However, their analysis involved 
enforcing the particular form of the spin-liquid order and 
imposing unequal coupling constants in Eqs. ^ and ([s]) 
to stabilize it (in our case they are both equal to Jh)- 




FIG. 4. KHLM at finite T. (a) Phase diagrams exhibiting 
competing orders, see text. Notations are the same as in Fig. 
[2] Shaded regions denote AF Kondo phases. The transition 
between Kondo and AF (AF Kondo) states is 1st (2nd) order. 
Unmarked phases are usual Fermi liquids, (b) OPs in the 
AF Kondo phase for a fixed Ju/t = 0.25 and t' /t = 0.9: 
(xo) (main panel and lower inset), and (xz) and staggered 
magnetizations (upper inset). 



We provided a framework where the fully self-consistent 
treatment of the KHLM allows controlled analysis to 
the different OPs. Crucially, while we find the non-zero 
"spin-liquid" order parameter in all the Kondo screened 
phases, it vanishes in the local moment antiferromagnet, 
and therefore the spinon FS is absent in that phase. 

It is known that the finite-T HMF phase transition be- 
tween HF and small-FS metal phases becomes a crossover 
when fluctuations beyond HMF are taken into account. 
However, the condensation of bosons xo and x does re- 
main a phase transition at T = [l2j, and hence salient 
features of the quantum phase diagram of Fig. [2ja) re- 
main unchanged, although phase transition lines may 
shift. At finite T we expect a wide quantum critical 
regime above those transition lines. 

We focused on commensurate magnetic phases. How- 
ever, at least in the classical KHLM such states are of- 
ten suppressed in favor of incommensurate spiral |21j or 
skyrmion | 22| phases. Therefore it would be interesting 
to establish whether such states are realized in the quan- 
tum KHLM. We leave this for future investigations. 
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